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I.  INTRODUCTION 


When  a  laser  beam  is  used  to  illuminate  a  rough  surface,  the  result¬ 
ing  scatter  pattern  takes  on  a  random,  mottled  character.  This  pattern 
is  usually  referred  to  as  a  speckle  pattern,  and  its  statistics  relate 
not  only  to  the  particular  radiance  pattern  on  the  surface  of  the  scat- 

terer  but  to  the  surface  roughness  as  well.  For  another  investigation1, 
a  computer  program  was  written  to  evaluate  numerically  the  far- field 
scatter  pattern  due  to  coherent  illumination  from  a  one-dimensional 
surface  characterized  by  gaussian  statistics.  For  that  particular 
study,  the  statistics  of  the  speckle  pattern  were  examined  with  respect 
to  the  form  of  the  radiance  distribution  over  the  scattering  surface. 

This  present  investigation  was  prompted  by  a  problem  related  to 
gun  erosion.  Central  to  the  estimation  of  gun  tube  wear  is  the  assess¬ 
ment  of  surface  roughness  from  one  firing  to  the  next.  With  a  view 
toward  developing  an  experimental  technique  useful  in  the  analysis  of 
surface  structure,  we  utilized  our  one-dimensional  analytical  model  to 
examine  the  relationship  between  roughness  at  a  scattering  surface  and 
certain  statistics  of  the  scattered  radiation. 


II.  PROBLEM  FORMULATION 

An  attempt  was  made  to  model  a  scattering  geometry  which  might  be 
used  in  an  actual  measurement.  With  this  in  mind,  an  object  1  mm  in 
length  was  chosen  with  radiation  at  0.6-micron  wavelength.  The  geometry 
is  shown  in  Figure  1.  The  line  source  formed  the  input  tc  a  15-cm 
focal  length  lens.  At  the  output  focal  line,  the  speckle  field  was 
evaluated  using  the  Huygens- Fresnel  equation.  The  integral  form  of 
this  equation  describes  the  output  field  as  an  integral  over  the  inpu\ 
aperture,  where 

X 


I, 

f - *4 

Figure  1.  Setup  for  generation  of  optical  scattering  pattern.  Laser 
light  impinges  on  surface  at  5  plane.  Radiation  is  transformed  by  lens 
to  plane  X.  The  problem  is  to  relate  the  statistics  of  the  optical 
radiation  at  plane  X  to  the  surface  characteristics  in  plane  c, 


1.  R.  Barakat,  P.H.  Deitz  and  T.E.  Buder,  "A  New  Class  of  Communica¬ 
tions  Systems:  II.  Theory",  BRL  Memorandum  Report  to  be  published. 
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V(x)  «f  V(^)expf-i  x  5]  d£,  (1) 

•'-a  «■  J 

Here,  V(x)  is  the  output  field,  the  limits  of  integration  are  from  -a 
to  a  where  a  is  h  mm,  V(£)  is  the  complex  scattered  field  at  the  rough 
surface,  X  is  the  wavelength  of  the  radiation,  and  f  is  the  focal  length 
of  the  lens. 

In  our  evaluation,  a  special  subroutine  was  used  which  generates 
stochastic  functions  which  are  gaussian  in  nature  and  conform  to  chosen 
parameters  of  standard  deviation  (vertical  surface  structure)  and  corre¬ 
lation  interval  (avei'age  scale  size  along  the  surface) .  Since  the 
illuminated  spot  is  small,  a  constant  amplitude  was  assumed  within  the 
scattering  aperture  so  that  the  statistics  of  the  scattered  radiation 
were  a  function  of  phase  only.  Figure  2  shows  one  computer  realization 
of  the  phase  of  the  optical  beam  immediately  after  scatter  at  the  rough 
surface  for  one  set  of  parameters.  The  standard  deviation  is  one-tenth 
of  a  wavelength  and  the  correlation  interval  is  5  microns. 


X  (mm) 


Figure  2,  Computer-generated  guasslan  phase  surface  having  standard 
deviation  (sigma)  of  0.1  lambda  (where  lambda  (l)  is  0.6  microns)  and 
a  correlation  interval  (alpha)  of  five  microns. 


Equation  (1)  was  utilized  in  discrete  form  for  400  equally  spaced 
points  across  the  1-mm  object  and  evaluated  at  40o  locations  in  the 
detection  plane  over  a  distance  of  40  cm.  After  the  400  complex  field 
calculations  were  made  in  the  detection  plane,  th<»  results  were  squared 
to  give  the  measurable  quantity,  the  intensity  pattern. 
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For  tho  gun  erosion  application,  a  range  of  height-roughness 
variations  from  0  to  0.9  X  (in  standard  deviation),  whore  X  is  the 
wavelength,  wore  examined  together  with  correlation  interval  variations 
from  5  to  10S  microns.  Six  particular  values  of  standard  deviation 
wore  used  with  five  values  for  tho  correlation  interval  in  every  possi¬ 
ble  pairwise  combination.  For  each  of  those  30  pairs  of  parameters, 

15  realizations  of  the  radiation  pattern  in  tho  detection  plane  wore 
generated.  Various  statistical  parameters  were  extracted  from  these 
speckle  patterns  and  are  discussed  in  tho  noxt  soction, 


II r,  RESULTS  OF  CALCULATIONS 

Using  the  15  realizations  of  speckle  pattern  for  each  of  30  pairs 
of  surface  height  (sigma)  correlation  interval  (alpha),  four  specific 
statistical  variables  wore  examined:  moan  irradianco,  standard  devia¬ 
tion,  ratio  of  RMS/nverago,  and  the  correlation  interval. 


For  each  of  the  four  statistical  parameters  examined,  the  pro¬ 
cedure  was  similar.  The  parameter  under  evaluation  was  computed  for 
each  speckle  pattorn  and  then  averaged  over  tho  15  realizations  cemented 
for  each  of  the  30  pairs  of  sigma/alpha  variables.  Standard  deviations 
about  those  moans  wore  also  computed  to  givo  an  indication  of  the  sta¬ 
tistical  uncertainty,  F.iguro  3  gives  the  oxpocued  value  of  tho  speckle 


Figure  3,  Three-dimensional  perspective  plot  showing  expected  value  of 
average  irradianca  of  scattering  pattern  vs.  pairwise  combinations  of 
surface  height  (sigma)  and  correlation  interval  (alpha)  of  illuminated 
surface.  L  is  the  wavelength  of  the  optical  Illumination  (0.6  microns). 
Each  point  represents  the  expected  value  of  15  realizations  of  a  speckle 
pattern  for  which  the  average  value  has  been  computed.  Error  bars  indi¬ 
cate  the  one-sigma  variations  about  the  expected  value.  Absence  of  error 
bars  Indicates  a  nonstochastic  calculation. 
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Figure  4.  Perspective  plot  showing  expected  value  of  standard  deviation 
vs,  pairwise  combinations  of  surface  height  (sigma)  and  correlation  inter 
val  (alpha)  of  scattering  surface. 


Figure  5.  Perspective  plot  showing  expected  value  of  RMS/AVE  ratio  v 

of  surface  height  (sigma)  and  correlation  inter 
(alpha)  ot  scattering  surface. 


Figure  6.  Perspective  plot  showing  expected  value  of  correlation  inter¬ 
val  (in  nra)  vs.  pairwise  combinations  of  surface  height  (sigma)  and  corre¬ 
lation  interval  (alpha)  of  scattering  surface. 

irrndinnco  (in  relative  units)  versus  various  sigma/ a. phas.  As  might  ho 
expected,  the  moan  is  essentially  independent  of  surf  ice  structure. 

Bars  above  and  below  the  mean  values  indicate  the  one-sigma  values  based 
on  IS  realizations.  An  increase  in  the  correlation  interval  gives  a 
greater  uncertainty  in  the  value  of  the  expocted  irradiance.  Absence 
of  error  bni's  indicates  a  nonstochastic  calculation. 

Figure  <1  gives  the  results  for  the  standard  deviation  of  speckle 
statistics  versus  the  pair-wise  independent  parameters ,  Standard 
deviations  are  generally  higher  for  low  surface  standard  deviation  and/or 
largo  correlation  interval.  Figure  5  shows  the  ratio  of  RMS  surface 
roughness  to  average.  This  ratio  follows  tho  same  trend  as  the  standard 
deviation  given  in  Figure  4.  This  result  is  to  be  expected  «unce  the 
standard  deviation  and  RMS  are  essentially  identical  computations  for 
a  large  sample  size,  and  the  mean  here  is  relatively  constant. 

Finally,  Figure  6  givos  the  expected  correlation  interval  tin  mm) 
as  a  function  of  the  sigma-alphas.  This  function  goes  through  a  notice¬ 
able  peak  for  a  smull  correlation  interval  of  surface  roughness  and  the 
midrange  of  surface  standard  deviation. 
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IV.  CONCLUSIONS 


The  surface  shapes  suggested  by  figures  4-6  indicate  that  it  may  be 
possible  to  associate  specific  scatter  statistics  with  particular  sur¬ 
face  parameters.  The  central  problem  is  to  associate  a  unique  set  of 
surface  conditions  with  a  particular  statistic  in  the  measured  speckle 
pattern.  From  these  results,  it  can  be  seen  that  specification  of, 
for  example,  a  particular  ratio  of  RMS/AVE  value  does  not  relate  to 
just  one,  but  a  restricted  set  of  alpha- sigma  pairs.  However,  the 
combined  use  of  both  the  RMS/AVE  results  (Figure  5)  together  with  the 
correlation  interval  calculations  (Figure  6)  would  help  reduce  the 
number  of  sigma/alpha  pairs  which  can  result  in  the  same  scattering 
statistics.  Further,  if  it  can  be  argued  on  physical  grounds  that  the 
correlation  interval  and  the  surface  height  of  the  scattering  surface 
musr  be  correlated  in  some  way,  then  the  ambiguity  in  the  interpreta¬ 
tion  of  the  scattering  data  might  further  be  reduced. 

From  such  optical  scattering  experiments  as  those  described  here, 
it  might  then  be  possible  to  relate  a  series  of  speckle  measurements 
to  a  specific  set  of  surface  roughness  characteristics. 
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150  Totwa  Road 
Wayne,  NJ  07470 
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DISTRIBUTION  LIST 


No.  of 

Copies  Qrgahanization 

1  System  Planning  Corporation 
Suite  1500 

ATTN:  Dr.  James  Meni 
1500  Wilson  Boulevard 
Arlington,  VA  22209 

1  University  of  California 
LawrCicc  Livermore  Laboratory 
ATTN:  Mr.  R.  Singleton 
P.0.  Box  808 
Livermore,  CA  94550 


Aberdeen  Proving  Ground 

Dir,  USAMSAA 

ATTN:  Dr.  J.  Sperrazza 

DRDXSY-MP,  Mr.  II.  Cohen 
9RXSY-GI,  Mr.  W.  Clifford 
DRXSY-DS,  Mr.  J.  Kramar 
DRXSY-T,  Mr.  A.  Reid 
DRXSY-0,  Mr.  R.  Conroy 
Mr.  C.  Odom 

Dir,  MEL 

ATTN:  Mr.  J.  Torre 

Mr.  G.  L.  I  lor  ley 
CDR,  USATEC0M 
A'lTN:  DRSTI5-  T0-I-* 

Dir,  Wpns  Sys  Concepts  Team 
Bldg  133516,  13A 
A’lTN :  DRDAR-AGW 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet  and  return  it  to  Director,  US  Army  Ballistic  Research 
Laboratory,  ARRADCOM,  ATTN:  DRDAR-TSB,  Aberdeen  Proving  Ground, 
Maryland  21005.  Your  comments  will  provide  us  with  information 
for  improving  future  reports. 

1.  BRL  Report  Number _ _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  rolated 
project,  or  ether  area  of  interost  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procodui'o,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  tho  information  in  this  report  led  to  any  quantitative 
savings  ns  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


S.  Gonoral  Commonts  (Indicate  what  you  think  should  be  changed  to 
mako  this  report  and  future  roports  of  this  typo  more  responsive 
to  your  needs,  more  usable,  improve  roadability,  etc.) _ 


6.  If  you  would  like  to  bo  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  tho  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 


